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I consider variations in the definition of a General-Mass Variable Flavour Number Scheme (GM- 
VFNS) for heavy flavour structure functions, both at next-to-leading order (NLO) and at next- 
to-next-to leading order (NNLO). I also define a new "optimal" scheme choice improving the 
smoothness of the transition from one flavour number to the next. At both NLO and NNLO I 
investigate the variation of the structure function for a fixed set of parton distribution functions 
(PDFs) and also the change in the distributions when a new MSTW-type global fit to data is 
performed for each GM-VFNS. At NLO the parton distributions, and predictions using them at 
hadron colliders, can vary by ^ 2% from the mean value. Use of the the Zero-Mass Variable 
Havour Number Scheme, which is simpler but only an approximation, leads to results a further 
couple of percent or more outside this range. At NNLO there is far more stability with varying 
GM-VFNS definition. Typical changes in PDFs and predictions are less than 1%, with most 
variation at very small x values. This demonstrates that mass-scheme variation is an additional and 
significant source of uncertainty when considering parton distributions, but like other theoretical 
uncertainties, it diminishes quickly as higher orders are included. 
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The tieatment of heavy flavours, charm and bottom, in structure functions has an important 
impact on the PDFs extracted in fits, due to direct data on F2{x,Q^) and also the contribution to 
the total structure function at small x. There are two distinct regimes with different descriptions. 
For massive quarks are created in the final state, and described using the Fixed Flavour 

Number Scheme (FFNS) (see [|l|] for NLO results), F{x,Q^) = cl^'"^ {Q^ /ml)^ fl^ {Q^), where «/ 
is the number of light quarks. This does not sum a^ln" Q^/ml terms in the perturbative expansion 
which may be important. At high scales, » m^, heavy quarks behave like massless partons and 
the logs are summed via evolution equations. The distributions for different light quark number 
are related to each other perturbatively = Ajk{Q^ /m^) ® f'^^ {Q^), where the matrix 

elements Ajk{Q^ /mjj), calculated at ff{ccj) in [Q], contain the fixed-order ln(2^/m^) contributions. 
In the Q^/ml — > oo limit the description is the Zero-Mass Variable Flavour Number Scheme (ZM- 
VFNS), F{x,Q^) — (j^^^'^'"/ ^fy{(2^). This is approximate, ignoring all &{m\/Q^) corrections. 
To correct this shortcoming and obtain a correct description between the two limits of < mjj 
and ^ mjj, one can use a General-Mass Variable Flavour Number Scheme (GM-VFNS). 

The GM-VFNS can be defined from equivalence of the ny flavour and + 1 flavour descrip- 
tions at all orders, resulting in d/'"' {Q^ /ml) = C™^'''"^^\q^ /ml) ®Ajk{Q^/ml), e.g. at ^(ttj) 

The VFNS coefficient functions tend to the massless limits as /ml oo, but Cj'^^ {Q^ /ml) 
is only uniquely defined in this limit. One can swap Gim^jQ^) terms between C^j^^''*'^(2^/m^) 

and C^^^^'^^^ in Eq. ([l]), and at higher orders, leading to various prescriptions [|3[, ^, 

^ III, 0]. The TR GM-VFNS highlighted the freedom in choice, and enforced correct kine- 
matics via a quite complicated definition. The (S)ACOT(;|^) prescription apphed the sim- 
ple choice C^j^^^''°^(eV'«I>z) °^ 5(z-Xmax), which gives f!1'^\x,Q^) oc el{h + h){x/x^^^,Q^), 
where x^ax = Q^/{Q^ and imposes the threshold = Q^{1 — x)/x> 4ml. gives 

the usual limit C^'^ '^'^'^\z) = 5{l—z) for /ml — > oo. The TR' scheme adopted this and 
extensions to higher orders (though uses a different multiplicative factor of Q^/{Q^+4ml) [p). 
However, ACOT-type schemes have used the same order of as above and below = ml, despite 
the fact that FFNS is LO at /^{as) while ZM-VFNS starts at zeroth order. Instead the TR' definition 
uses, for example, at LO the &{(Xs) FFNS result for < ml, and for > m| 

F,\x,Q') = as{ml)cZ:''^'\l)^ff{ml) +C^^'^"''^'^^'\QVml) {h + m'l (2) 

i.e. it freezes the higher order as term when going upwards through = ml. This difference in 
choice can be phenomenologically important. As an alternative, but ultimately equivalent formu- 
lation of a GM-VFNS, BMSN ^ and FONLL [§] define a scheme in general terms as 

where the second (subtraction) term is the asymptotic version of the first, i.e. all terms i^{ml/Q^) 
are omitted. There are differences in exactly how the second and third terms are defined in detail 
in different schemes. In the simplest applications the as order of F^{x, Q^) at low is the same 
as that of F™^^{x, Q^) as oo. There is a version of FONLL which uses one power higher in 
the FFNS term, but it leads to part of the higher order contribution persisting as — > oo. 
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Figure 1: The variation in F2 (x, 2^) generated from a variety of choices of GM-VFNS at NLO (left) and 
NNLO (right) using the MSTW2008 pdfs in each case. 

Ideally one would like any GM-VFNS to reduce to 
exactly the correct order FFNS at low and exactly 
the correct order (one power of as lower) ZM-VFNS 
as — ^ 00. At present none do, but this can easily be 
rectified. Let us return to the TR' version of the GM- 
VFNS. The obstacle is the presence of the frozen term 
as ^ 00 (which depends on the PDFs only at low 
scales, so is a small effect at large Q^). In fact, this is 
not strictly necessary and one can have instead 
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Table 1: Parameter values for different ex- 
ti-eme GM-VFNS definitions. 

(ml/Q^ra^imDl^C^Mmi'nf,) or (ml/Q^y a'iiQ^)l^C^,iQ')mQ')- (4) 

Any a > provides the correct limit, though strictly from factorization one should have {mjj/Q^) 
times \n{Q^/mjj) terms. There is also more freedom. One can modify the heavy quark coefficient 
function as long as the Q^/ml^°° limit is maintained. However, since this appears in convolutions 
for higher order subtraction terms, we do not want a complicated x dependence. A simple choice is 



c 



GMVF,{0) 



{Q'/mi,z) ^ (l+^K/G")^)5(z-x„,ax) 



(5) 



where again variation in c really mimics {m\/Q^) with logarithmic corrections. One can also 
modify the argument of the 5-function, similar to the Intermediate-Mass IM scheme [§], 



■x{\ + {x{\+Aml/Q^))Hml/Q^) 



(6) 



so the kinematic limit stays the same, but if J > 0(< 0) small x is less (more) suppressed. The 
default a,b,c,d are all zero, but can vary, being limited by fit quality or sensible choices. 

A variety of different choices defined in Table [T] has been tried at NLO and at NNLO along 
with the ZM-VFNS (at NLO). The resulting variations in F2{x, Q^) near the transition point due to 
different choices of GM-VFNS at NLO are shown in the left of Fig. [T[ I also define an "optimal" 
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Figure 2: The variation of PDFs obtained from the best fit from a variety of choices of GM-VFNS and the 
ZM-VFNS at NLO (left) and NNLO (right) as a ratio to the MSTW2008 PDFs. 

scheme which is smooth at threshold and reduces to exactly the right limits at high and low Q^. 
There is quite a spread at NLO, though the ZM-VFNS is far steeper at low than any GM-VFNS. 
This spread is very much reduced at NNLO, the right of Fig. |l|, with almost zero variation until 
very small x, showing that NNLO evolution effects are most important in this regime. 

Global fits are also performed using the 
same procedure as the MSTW08 fit [0] for 
all schemes. At NLO the initial for a new 
scheme can change by 250, but converges 
to within 20 of the original. There are im- 
proved fits for options 1, 3 and 6 and the 
fit is best for the for optimal scheme. The 
variations in the partons extracted at NLO 
are shown in the left of Fig. ^ The default 
TR' scheme sits near the low end. Some 
changes in PDFs exceed the one a uncer- 
tainty. a5(M|) changes by < 0.0007 except 
for the ZM-VFNS where it falls by 0.0015. 
The ZM-VFNS PDF is clearly outside than 
the GM-VFNS band. For fits at NNLO the 
initial changes in are < 20 and they converge to within 10 of the original. The variations in 
PDFs extracted at NNLO are shown in the right of Fig. At worst the changes approach the 
uncertainty, but are usually far less. Variations in as{M^) are ~ 0.0003. However, at NNLO the 
TR' scheme models the ^ _ _ 

)'(ln(l/z) 
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+ 1.1% 
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-0.1% 
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-0.3% 
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+ 1.6% 


+0.8% 


GMvaropt 


+0.3% 


-1.5% 


+2.0% 
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ZM-VFNS 


-0.7% 


-1.2% 


-3.0% 


-3.1% 


GMvarcc 


+0.0% 


-0.1% 


+0.0% 


-0.1% 



Table 2: Predicted cross-sections at NLO for Z and a 
120 GeV Higgs boson at the Tevatron and LHC. 



(a^) FFNS terms at low Q using leading threshold logarithms [11] and 



-b)/z, where the default is 



ln(l/;c) terms 1 12]. The latter take the form oc (1 
a = 20, ^ = 4, so a and b can be varied. Changes in a make little difference. The maximum sensible 
variation to ^ = 2 leads to an effect of order the uncertainty at .x < 0.001. However, this is largely 
eliminated if the ^(o;|) contribution dies away, rather than being frozen. 
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The predictions for cross-sections aie 
shown at NLO in Table. ^ There is at most 
a 1.5% variation at the Tevatron. There is 
a +3% down to —0.5% variation in Oz at 
the LHC. The spread in Oh is about halved 
due to the higher average x sampled. The 
ZM-VFNS is the clear outher in the low di- 
rection at the LHC. GMvarcc denotes vari- 
ation in the GM-VFNS for charged current 
processes, and clearly the effect is very 
small indeed. The predictions at NNLO 
are seen in Table. ||. Other than model de- 
pendence - GMvarmod denotes the varia- 
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tion to ^ = 2 in the ^{oc^) term - the max- Table 3: Predicted cross-sections at NNLO for Z and a 
imum variations are of order 0.5% at LHC. 120 GeV Higgs boson at the Tevatron and LHC. 
GMvarmod' is when the ^(a|) terms fall with Q^, and also exhibits a very small deviation. 
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